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Abstract. The precision measurement of the 1809 keV gamma-ray line from Galactic 26 Al is one of the goals of 
the SPI spectrometer on INTEGRAL with its Ge detector camera. We aim for determination of the detailed shape 
of this gamma-ray line, and its variation for different source regions along the plane of the Galaxy. Data from the 
first part of the core program observations of the first mission year have been inspected. A clear detection of the 
26 Al line at ~ 5-7 a significance demonstrates that SPI will deepen 26 Al studies. The line intensity is consistent 
with expectations from previous experiments, and the line appears narrower than the 5.4 keV FWHM reported 
by GRIS, more consistent with RHESSI's recent value. Only preliminary statements can be made at this time, 
however, due to the multi-component background underlying the signal at ~ 40 times higher intensity than the 
signal from Galactic 26 Al. 
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1. Introduction 



The measurement of 1809 keV emission from Galactic 
26 Al has been one of the design goals of the INTEGRAL 
mission l|Winkler et al.U2003tlHermsen fe Winklerl.l2002jh 
26 Al gam ma-rays were discove red already in 1982 by 
HEAO-C l|Mahonev et al.Lll982f) . and since then are con- 
sidered direct proof of ongoing nucleosynthesis in the 
Galaxy. Several follow-up experiments have set out to 
measure details about 26 Al sources. 

Much has been learned through the Compton Gamma- 
Ray Observatory, in particular with COMPTEL's sky 
survey over 9 years, which resulted in an all-sky im- 
age in the 1809 keV gamma-ray line llPluschke et al 



200lUKnodlseder et allfl999tlOberlacklll997t Diehl et al 
1995|b This image clearly demonstrates that emission 
along the plane of the Galaxy dominates, hence 26 Al nu- 
cleosynthesis is common throughout the Galaxy, rather 
than a local phenomenon of the solar system. The irreg- 
ular structure of the emission and alignments of emission 
maxima with spiral-arm tangents sugg ested that massive 
stars dominate 26 Al nucleosynthesis ijChen et all Il995l 
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iPrantzos fe DiehlL I 1996). This could be further substan- 
tiated through comparisons with candidate source trac- 
ers, and through modelling of 26 Al emission from the 
Galaxy and specific source regions ba sed on knowledge 
about the massive-sta r populations ijKnodlseder et all 
ll999llKnodlsederMl999lh 

The high spectral resolution of Ge detectors of 3 keV 
(FWHM) at the 26 Al line energy of 1808.7 keV is ex- 
pected to reveal more information about the sources and 
their location through Doppler broadenings a nd shifts, 
from Galactic rotation l|Gehrels fe ChenL Il996li and from 
dynamics of the 26 Al gas ejected into the interstellar 
medium. In particular after the GRIS balloon experi- 
ment and their re port of a significantly-broadened line 
( Nava et all fl996i. alternative measurements of the 26 Al 
line shape were of great interest. GRIS's value trans- 
lates into an intrinsic line width of 5.4 keV, equivalent 
to a Doppler broadening of 540 km s _1 . Considering the 
1.04 x 10 6 y decay ti me of 26 Al such a l arge line width is 
hard to understand ijChen et all Il997t ISturner fe Naval 
1999). Other high-resolution measurements are in mild 
conflict with the GRIS result. The original HEAO- 
C measurement was interpreted as an intrinsically- 
narro w line of width less than 2 keV ijMahonev et all 
Il984h . and the recent RHESSI measurement shows some 
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broadening, however at an inter mediate value of about 
2.0(4-0.78,-1.21) keV FWHM l)Smitbl . |2003|) . signifi- 
cantly less than the GRIS value. SPI on INTEGRAL 
l|Vedrenne et all l2003t fRoaues et all 12003^) with its com- 
petitive spectral r esolution and the INT EGRAL core ob- 
serving program ilWinkler et all 1200.1) emphasizing ex- 
posures of the inner region of the Galaxy is expected to 
clarify these questions through high-quality data. The ini- 
tial calibration phase early in the mission, which used 
exposures of sources in the Cygnus region, had demon- 
strated an excellent performance of the instrument and 
even shown convincing evidence of detections of diffuse 
26 Al emission from the Cygnus region, supporting these 
prospects. In this paper, we report initial analyses of 
the first part of the inner Galaxy region survey of the 
INTEGRAL core program. 

GCDE(l) Pointings Rev 47-66 
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Fig. 1. Pointings of the first part of the Galactic-Central- 
Deep-Exposure analyzed here 



orbit. Calibrated single-detector events ("SE", triggering 
one single of the 19 Ge detectors), and the composite of 
single and multiple detector hits (SE+ "multiple events, 
ME") were analyzed. Note that at the 26 Al line energy 
40% of the measured events are multiples. A systematic 
uncertainty of our energy calibration at energies below 
~ 150 keV may however distort the energies assigned to 
multiple events. We therefore perform independent analy- 
ses per event type (SE, SE+ME), to check for such system- 
atics which could lead to artifacts in the wings of spectral 
lines. 
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Fig. 2. Raw spectrum of Crab reference observations 
(SE). Background dominates the signal, with prominent 
instrumental lines at 1779 keV originating from 28 Al, at 
1764 keV from 205 Bi, and at 1809-1811 keV from a blend 
of 26 / 27 Na and 56 Mn. 



2. Observations, Data, and Analysis 

The first year of the INTEGRAL mission will emphasize 
a deep survey of the inner part of the Galaxy, devoting 
about 4 Msec of exposure to this region, in a dedicated 
observing pattern characterized by a 2° pitch, extend- 
ing ± 30° in longitud e and ± 20° in latitu de away from 
the Galactic Center (I Winkler et all [2003). Data shar- 
ing agreements within the INTEGRAL Science Working 
Team imply that results on Galactic 26 Al for the sepa- 
rate Galactic quadrants will not be addressed here, but 
reported in future papers. The data used for this ini- 
tial analysis encompass INTEGRAL orbits 46-66 (mis- 
sion days 1157-1216). Not all our analyses use the exactly 
same data, different selections were applied from early 
stages of processing; a minimum set comprises in total 
839 pointings with 0.971 Ms integrated exposure livetimc 
(see Figure 0. 

Energy calibration during this time was derived by fit- 
ting instrumental background lines at energies 439, 585, 
883, 1014, 1369, and 1779 keV as accumulated for each 
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Fig. 3. Raw spectrum of Crab reference observations for 
multiple events (ME). The same instrumental lines are 
seen, slightly weaker if compared to continuum back- 
ground. 
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Fig. 4. Raw spectrum of GCDE measurements (SE). The 
form of the instrumental-background spectrum is virtually 
identical to the Crab reference (Fig.|2Jl, the signal expected 
from celestial 26 Al is at the percent level. 



Background dominates the overall signal, (see Figures 
EffiJ, so that in the 26 Al line region (1809±4 keV) we mea- 
sure 5 x 10 5 counts from the inner Galaxy, while our back- 
ground reference is based on exposures of the Crab region 
providing 1.3 x 10 5 counts in the line-region. Continuum 
dominates, but ~ 17% of the total signal is in a rather 
broad line-like feature around 1810 keV. This is mostly 
instrumental background, which needs to be understood 
before extracting the 26 Al signal: For the GCDE, we ex- 
pect from COMPTEL measurements a signal strength 
of about 13000 counts or 2.7% of the total measured 
counts in the line region. Background is expected at 1808.7 
keV from excited 26 Mg produced from spallation of Al 
and from a-captures on Na, and at 1810.7 keV from 
56 Mn(/3 ~) 56 Co(.EC) 56 Fe, but other nu clear lines may con- 
tribute l|WeidensDointner et all [2003(1 . Correlation anal- 
ysis with other line features and cosmic-ray activation 
monitors is underway to model details of this underlying 
backg round, both in shape and intensity fe.g. lJean et all 
120031) . 

Data analysis is complicated due to the large number 
(839) of individual pointings, different measurement times 
of source and background, and detector and background 
evolution within and between them. For a strong source, 
one may hope to subtract a sufficiently well-defined back- 
ground and then sec a source signal. Indeed, when we use 
exposures from the Crab region, adjust for the different ex- 
posures and detector resolutions by normalization on the 
nearby background lines at energies 1764 and 1779 keV, 
and subtract this normalized "off-source" reference from 
the integrated spectrum measured from the inner Galaxy, 
we obtain an excess signal which indicates the presence 
of a celestial signal from 26 Al (Fig[5j). The strong instru- 
mental background lines do not perfectly subtract. The 
suppression by a factor ~ 90, however, would correspond 
to residuals from the instrumental feature at 1809 keV 



of less than half of what we observe; also, residuals ap- 
pear in the wings of the instrumental lines, because our 
normalization does not trace time-variable gains or non- 
Gaussian detector degradation. The feature at 1809 keV, 
however, resembles the expected line more closely than 
background residuals: we indicate the expected 26 Al signal 
with a Gaussian at the 26 Al line energy and instrumental 
line width. Note that here we do not use SPI's imaging ca- 
pabilities, so it is not surprising that a not very significant 
26 Al detection is obtained (see Figure EJ). 

Imaging analysis makes use of the detailed response 
of the SPI i nstrument as is obta ined from Monte Carlo 
simulations l(Sturner et all 120031) . including the proper- 
ties of the coded mask, the Ge camera, and all other 
material possibly interacting with incident gamma-rays . 
Calib rations before INTEGRAL'S launch jAttie et all 
2003J) and on Cyg-Xl and the Crab in the first part of 
the mission have verified the validity of this response 
l|R.oaues et all 120031 ISturner et all l2003h . Adopting an 
intensity distribution on the sky, one may thus predict the 
measured event distribution from 26 Al for the observed 
exposure pattern. Adding these to a suitable background 
model, one obtains a prediction for the actually-measured 
dataset. A fitting procedure can then be used to adjust 
intensity parameters of the 26 Al skymap and background 
components: 



D 



i,3,E 



Ie ■ S ® 3? + a hE ,k ■ B E ,k + R ijiE 



Here D are the measured data, S the sky intensity distri- 
bution, 5R the instrument response function, B the back- 
ground model, H the statistical noise, and indices are i for 
detectors, j for pointings, E for energy, and k for back- 
ground model components. From such fits we obtain in- 
tensity spectra Ie of our sky model as constrained by 
our measurement. Several implementations (e.g. spidif- 
fit/spi obs-fit) of thi s analysis concept have been devel- 
oped ijStrongl l2003t IKnodlsederl l2003h . differing in the 
method of minimum searches and uncertainty estimates: 
The fit determines the mean of the intensity parameter 
posteriors through a Markov-Chain Monte Carlo method 
marginalizing over the background parameters (spidif- 
fit), or alternatively minimizes the log- likelihood function 
through a Levenberg-Marquardt minimum search to fit 
spectral intensities (spi_obs_fi,t) . Background amplitudes 
o-j.E.k are fitted per pointing and energy from the measure- 
ment (spidiffit), or prescribed by an absolute model based 
on tracers of continuum and line background components 
(spLobs-fit) . Error bars SIe are determined as standard 
deviation from the mean with the same Markov-Chain 
Monte Carlo method fetrong L l2003|) for the large number 
of fitted parameters in spidif fit, and with an eigenvalue 
analy sis of the error matrix <|StrongL Il985t IKnodlsederl 
2003) for the few fitted parameters in spi_obs-fit. 

Estimates of the systematic uncertainties have been 
added in quadrature to these statistical uncertainties and 
are thus included in the error bars given in figures 
IH1 We fit different models (fc) of background, such as 
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continuum background interpolated from adjacent ener- 
gies, line backgrounds from scaled reference observations 
off the source of interest, and scaling models using suit- 
able tracers of background such as the rate of saturated- 
signal events in the Ge detectors; these are assumed 
to arise from cosmic-ray triggers which activate space- 
craft material and hence generate nuclear-line background 
l|jean et all 12003^ . Suitable models S of the 26 Al sky in- 
tensity distribution were adopted f rom t he COMPTEL 
26 Al all-sky results l|Pluschke et alll200l|) . and from dis- 
tributions of free electron s or warm dust as derive d from 
COBE measurements (see iKnodlseder et all Il999h . 

If the signal is sufficient, one may consider imag- 
ing analyses in narrow spectral bins, such as exercised 
for the case of the 511 keV e mission from the inner 
Galaxy l|Knodlseder et all l2003(k in our case, the signal 
is too weak to reasonably constrain such analysis with 
its intrinsically many free parameters. For instrumental 
background only, no or only modest correlation with a 
plausible 26 Al skymap is expected, while the correlation 
should be significantly better within the 26 Al line energy 
bins. However, systematic uncertainties in our background 
models lead to distortions of such a correlation. It is the 
main task of current analysis efforts to understand the 
level and possible spectral shapes of such systematic un- 
certainties. 
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Fig. 5. Residual spectrum of GCDE - Crab (SE), after 
normalization on the intensities of the 1764 and 1779 
keV instrumental background lines, adjusting for detec- 
tor degradation between the two observations. The resid- 
ual excess around 1809 keV appears above the properly- 
scaled residuals from the instrumental background lines, 
which remain from imperfections in the subtraction pro- 
cess. Here the coded-mask imaging information is not 
used. 



line, we obtain clear detections of celestial 26 Al emission 
at the level of 5-7cr. The results for the 26 Al flux, as well 
as details of the spectral signature, however, vary signifi- 
cantly with parameters of the analysis, and thus indicate 
the levels of uncertainty at this initial stage of the work; 
statistical uncertainties are negligable, in comparison. 
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Fig. 6. Imaging analysis result from fitting a sky inten- 
sity distribution as modelled from the COMPTEL 26 Al 
skymap to each energy bin. Background was modelled 
from Crab observation detector ratios, and fitted in inten- 
sity to the actual measurement for each pointing, together 
with the sky signal. 
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Fig. 7. Imaging analysis result from fitting a sky inten- 
sity distribution as modelled from the COBE/DIRBE dust 
skymap (240/im) to each energy bin. Background was 
modelled from adjacent energies for continuum, the Crab 
observation detector ratios were used for the line, and 
scaled by the rate variations of saturated events in the 
detectors to model the line background. 



3. Results 

From spectral analysis through fitting of adopted models 
for the 26 Al skymap over an energy range around the 26 Al 



In Figure|H|we show a spectrum derived from all (single 
and multiple) events with spidiffit using the COMPTEL 
Maximum-Entropy map from 9 years of measurements as 
a model for the spatial distribution of the sky emission 
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Fig. 8. Same as Figure but for the exposures of the 
LMC/SN1987A used as background reference. 





FWHM [keV] 


I [lCT 4 ph cm" 2 s _1 l 


uncertainty 


0.7 


1.4 


fit value Fig. El 


2.1 


3.3 


fit value Fig. 


3.1 


3.3 


fit value Fig. |51 


3.1 


4.7 



Table 1. The range of values fitted for several analysis 
scenarios illustrates the range of systematic uncertainties. 



COMPTEL had measured a value somewhat lower value 
of 2.8 ± 0.15 x 1(T 4 ph cm^s" 1 l|Oberlackl Il997|k with 
their background subtraction from high-latitude observa- 
tions possibly suppressing large-scale diffuse flux compo- 
nents. On the other hand, RHESSI recently mea sured a 
rathe r high value of 5.7 ±0.54 x 10~ 4 ph cm _2 s~ 1 ijSmithl 
l2003h from earth occultation analysis of their measure- 
ments pointed at the sun. 

The different values for our line positions, all some- 
what higher than the expected value of 1808.7 keV, may 
indicate systematics in our energy calibration, the impact 
of detector degradation over the time of measurements, 
or some structured underlying background effect which 
may shift our signal upward by a few tenths of a keV. 
On the other hand, all line width results which we ob- 
tain are consistent with SPFs instrumental resolution of 
3 ke y (FWHM) a nd thus support RHESSI's recent find- 
ing ilSmithl.l2003h that the broad line reported by ORIS 
l)Nava et all Il99fih probably cannot be confirmed. Work 
is in progress to refine our spectral treatment and back- 
ground modelling, in order to be able to further quantify 
and substantiate this conclusion. 
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ijPluschke et all l200lh . Given the rather modest spatial 
resolution of SPI, the particular choice of such distribution 
is not critical, as long as the dynamic range of fluxes and 
spatial distribution are approximately correct; any choice 
of good source tracers, such as the warm dust or fre e elec- 
tron distributions (see e.g. iKnodlseder et all 119991) . pro- 
duce very similar results. Here we treat background by 
adopting the relative detector intensity ratios from the 
Crab exposures, and adjust their intensity for each point- 
ing in the fit. 

In FigureEJwe compare this to a spectrum derived from 
single events only (to avoid ME energy calibration issues), 
with spi.obs.fit and the COBE/DIRBE 240 ^m dust map 
as sky model, modelling background for the continuum 
from adjacent energy bins, and for the line component 
through taking detector ratios from the Crab exposures 
and absolute intensity variations from the rate of satu- 
rated events in the Ge camera during the actual GCDE 
measurements . 

Figure [3] then illustrates how results depend on differ- 
ent datasets for the background: here the exposures taken 
for SN1987A were used as a reference to model detector 
ratios, in an otherwise analogous analysis with spi.obs.fit. 

Our fitted sky intensity values (see Table 1) from the 
inner ±30° of the Galaxy are (3-5) xl0~ 4 ph cm _2 s _1 , 
and thus fall into the range suggested by previous observa- 
tions: The general consensus for the inner-Galaxy bright- 
ness of 26 Al is ~ 4 x 10~ 4 ph cm~ 2 s~ 1 , cons idering all 
uncertainties involved l|Prantzos fe Diehl I 1996); here "in- 
ner Galaxy" means integrating over the central radian, 
roughly ±30° in longitude about the Galactic Center. 
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